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Abstract 

We consider a heavy piston in an infinite cylinder surrounded by ideal gases on 
both sides. The piston moves under elastic collisions with gas atoms. We assume 
here that the gases always exert equal pressures on the piston, hence the piston 
remains at the so called mechanical equilibrium. However, the temperatures and 
densities of the gases may differ across the piston. In that case some earlier studies 
by Gruber, Piasecki and others reveal a very slow motion (drift) of the piston in 
the direction of the hotter gas. At the same time the hotter gas slowly transfers 
its energy (heat) across the piston to the cooler gas. While the previous studies of 
this interesting phenomenon were only heuristic or experimental, we provide first 
rigorous proofs assuming that the velocity distribution of the ideal gas satisfies a 
certain "cutoff" condition. 

1 Introduction 

Consider an isolated cylinder filled with an ideal gas and divided into two compartments 
by a large piston which is free to move along the axis of the cylinder, Fig. 1. The piston 
interacts with the gas atoms via elastic collisions. Assume that the gas in each com- 
partment separately is at equilibrium with temperature and density T_,n_ and T + ,n+, 
respectively. Let the gases exert equal pressures on the piston, i.e. let 

P_ = n-ksT- = n + kBT + = P + 

Then the system is at the so called mechanical equilibrium, and according to the laws of 
thermodynamics this state should be (macroscopically) stable. 
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Figure 1: Piston in a cylinder filled with gas. 



However, the system as a whole is not in a true equilibrium state, unless T_ = T+, 
hence microscopically it is not stable yet and should find ways to evolve to a true, thermal 
equilibrium, in which T_ = T+. This was predicted earlier by Landau and Lifshitz 
[LL], Feynman [F] and others. Recently Gruber, Piasecki and Frachebourg and others 
[GP, GF, P] derived heuristically, by means of kinetic theory, exact formulas describing 
the slow drift of the piston and the slow heat transfer from the hotter gas to the cooler 
gas. 

Our goal is to derive rigorously the main formulas of [GP, GF, P] describing the slow 
drift of the piston and the slow heat transfer between the gases. 

The piston model trivially reduces to a one-dimensional system by the projection onto 
the axis of the cylinder. Then one obtains an ideal gas on an interval, and the piston 
itself becomes a heavy point particle. The motion of a heavy particle in an infinite ideal 
gas of light particles is a classical example of Brownian motion studied by van Kampen 
[vK] and many others [L, H, DGL, GF, GP, P]. 

We consider a one-dimensional ideal gas on the entire line, without boundaries. The 
heavy piston is initially placed at the origin X(0) = and is at rest V(0) = 0. The initial 
configuration of gas atoms and their velocities is chosen at random as a realization of a 
(two-dimensional) Poisson process on the (x, t>)-plane with density p(x,v). This means 
that for any domain D C R 2 the number Np of gas particles (x, v ) & D at time t = is 
a Poisson random variable with parameter 



The system evolves according to the rules of elastic collisions. Denote the mass of the 
piston by M and the mass of an atom by m. Since atoms have identical masses, their 
collisions can be ignored. When an atom with velocity v collides with the piston, whose 
velocity is V, their velocities after the collision are given by 
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These rules preserve the total kinetic energy and the total momentum. Between collisions, 
all the particles and the piston move with constant velocities. 

The position X(t) and velocity V(t) = dX(t)/dt of the piston make a random process 
whose characteristics are determined by the initial gas density p(x,v). It is natural to 
assume that p(x, v) is symmetric in v and spatially homogeneous, i.e. p(x, v) = p(v) and 
p(v) = p(—v). We can set p(v) = nf(v), where f(v) = f(—v) is a probability density 
and n > is a (constant) spatial density. Then one can approximate X(t) and V(t) by 
certain Gaussian stochastic processes: 

Theorem 1.1 (Holley [H]) Let the density f(v) have a finite fourth moment j v 4 f(v) dv < 
oo. Then for every finite t < oo, the function V(t)y/~M on the interval [0, to] converges, 
in distribution, as M,n — > oo and M/n — > const, to an Ornstein-Uhlenbeck velocity 
process V t , while X(t)y/~M converges to an Ornstein-Uhlenbeck position process X t . 

An Ornstein-Uhlenbeck process (X t , V t ) is defined by [Ne] 

dX t = V t dt, dV t = -aV t dt + ^/D dW t 

where a > 0, D > are constants and Wt a Wiener process. The Ornstein-Uhlenbeck 
position process X t converges in an appropriate limit (e.g. a — > oo, a 2 / D = const) to a 
Wiener process. 

We note that Diirr et al [DGL] extended the above theorem to arbitrary dimensions. 

Our paper concerns with another physically interesting situation, where the initial 
density is not spatially homogeneous, but a spatial homogeneity is assumed separately 
for the gases to the right and to the left of the piston. So we assume that 

, s _ J p+(v) for x > 
P{x,v) - | p_( v j f or x < o 

and we also assume symmetry: P-(v) = p-(—v) and p+(v) = p + (—v). We can set 
P±( v ) = n ±f±( v ), where f±(v) = f±(—v) are probability densities and n± > are 
(constant) spatial densities. 

In addition, we want to exclude macroscopic motion of the piston in either direction. 
The macroscopic velocity of the piston vanishes if and only if 

ri- j v 2 f-(v)dv = n + / v 2 f + (v)dv (1.3) 
Jo Jo 

as it was shown heuristically in [LPS] and under some conditions rigorously in [CLS]. 
The physical interpretation of Eq. (1.3) is the pressure balance. If we define the pressures 
of the gases by 

/oo 
v 2 f±(v) dv 
-oo 

(which would be a proper thermodynamical pressure if the velocity distributions were 
Maxwellian) , then the condition (1.3) means exactly that P_ = P + . 
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To simplify some technical considerations we assume that the initial densities of the 
gases satisfy the following velocity cutoff: 

p±(v) = 0, if \v\ < v min or \v\ > v max (1-4) 

for some < i> m i n < f max < oo. Hence, the initial velocities of atoms are bounded away 
from zero and infinity. Under these conditions, our arguments are rigorous. We also 
discuss in Section 4 how to relax these assumptions, leaving this work for the future. 



2 Markov approximation 

Our cutoff condition (1.4) has an important implication. As long as the speed of the 
piston remains small enough, every gas atom collides with the piston at most once. 
Indeed, let \V(t)\ < V max , where 

M -m 

"max — „ , 7 ^min 

3M + m 

Note that V max is close to i> m in/3 when M ^> m. Then it follows directly from (1.4) and 
(1.2) that the atoms' velocities after collisions are at least 

M -m 2M 

" ^min „ r ' "max "max 



M + m M + m 

the last equality holds due to our choice of V max . Therefore, the atoms after collisions 
remain faster than the piston, so the latter cannot "catch up" with them. 

Hence, as long as |V(£)| < V ma , x , the velocity of the piston V(t) evolves as a Markov 
process with piecewise constant trajectories (a jump or step process). Moreover, this is a 
stationary (homogeneous in time) Markov process due to our requirements on the density 
p(x,v). We will slightly change the function V(t) so that it will evolve as a stationary 
Markov process unconditionally. Fix some V G (0, V max ) and require that whenever the 
velocity V of the piston after a collision, in the notation of (1.1), exceeds V in absolute 
value, it gets "reflected" at V, i.e. it instantaneously changes from V to V" by the rules 

V'>+V =}► V" = +2V-V ( 
V'<-V =^ V" = -2V-V [ ' 

These rules are in the spirit of random walks with reflecting boundary conditions. 

We will denote the velocity of the piston in the so defined dynamics by Wit). It is 
clear that Wit) = V(t) for alH < T such that sup 0<t<T \V(t)\ < V. We first study the 
Markov process Wit) and later estimate the difference V(t) — Wit). 

Denote by V(u, dw; At) for every At > the transition probability for the process 
Wit), i.e. 

P(W(t + At) G A/Wit) = u ) = [ V(u, dw; At) 

J A 
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for every Borel set AcR. It is clear that the piston does not experience any collisions 
with the atoms during the interval (t, t + At) if and only if the trapezoidal domain 

{v — u \ \ 

^'^ : X(t) < ~A? Vmin < ' U ' < VmBK \ ^ 

does not contain gas atoms. Therefore, the probability that W(t + At) = W(t) — u is 

V(u, {u}; At) = exp (- p(x, v) dx dv^j (2.3) 

To evaluate this integral, we partition the domain D = D(t,At) into two trapezoids, 
D = D~ U D + as shown on Fig. 2. Under our assumptions (1.4) 

f I" /""max r-Vmin 

/ p(x : v) dx — At / (v — u)p-(v) dv + / (u — v)p + (v) dv 

JD L^Vmin J~ Umax 

We introduce the following notation: for each k > let 

/•"max _ f-fmin , 

/ v p-(v) dv = F k and / v p + {v) dv = F£ 

t '' l? min *^ "^rriax 

and 

<2fc = / 'a- - F k 

Then we obtain 

/ p(x, v) dx = At \f~ - F+ - (F- - F+)u 

= At [Qi - Q u] (2.4) 

where u = W(t). 

It is clear that for every u and At > the probability measure V(u, dw; At), besides 
having an atom at u with probability given by (2.3), has an absolutely continuous com- 
ponent with a positive density on the interval (— V, V) (in fact, that density is bounded 
away from zero for every fixed u and At > 0). 

Proposition 2.1 The stationary Markov process W(t) has a unique stationary measure 
jiQ, which is absolutely continuous on the interval (—V, V). Any other initial distribution 
converges to p, uniformly exponentially fast in time. 

Proof. The process W(t) satisfies the so called Doeblin condition [D]: there exist e > 
and At > such that for every Borel set A C (— V, V) and every u G (— V, V) we have 



m 



(A) < e =^ f V(u, dw; At) < 1 - e 

J A 



where m(A) is the Lebesgue measure of A. Then our result follows from general theorems 
[D]. But we also outline a simple direct argument. 
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Figure 2: Region D = D~ U £>+. 

Since the transition probability T^w, g?u>; At) defines, for each At > 0, a continuous 
map on the convex space of probability distributions on the interval (— V, V), the exis- 
tence follows from a general Schauder-Tychonoff theorem [DS], p. 456. Since V(u, dw; At) 
has an absolutely continuous component whose density is bounded away from zero, the 
stationary distribution fi is absolutely continuous with density fo(w) > c > 0. 

To prove the uniqueness of /i , suppose /ij, ^ /i is another stationary measure with 
density f'^iw) > c' > 0. Then (1 + c)fo — cf for small c > is also a stationary density. 
Let 

c = sup{c > : inf_ [(1 + c)f (w) - cf' Q (w)\ > 0} 

|ui|<V 

Then (1 + c)f — cf is a stationary density, hence it must also be bounded away from 
zero on (— V, V). But this contradicts our choice of c. 

To prove the convergence, we take an arbitrary initial distribution z/ and consider its 
image u t at time t. Fix a t > 0. The measure u t has an absolutely continuous component 
whose density is bounded away from zero. Hence v t — (1 — a)u^ + ct/io with some small 
a > and some measure v^ 1 ' . The value of a depends on t but not on uq. This implies, 
by induction, that fkt = (1 — a) k u^ + [1 — (1 — a) k ]jjQ for all k > 1, which proves the 
exponential convergence of v t to /io- □ 

We now derive one useful equation. Fix a At > and for every u G (—V,V) denote 

by 

E(u,At) = I wV(u,dw;At) 
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the conditional expectation of the process W(t+At) given that W(t) = u. The invariance 
of the measure /i implies 



E(u, At) d/j,o(u) = / wdfj,o(w) 



Now suppose that 
so that 



E(u, At) = u + g(u) At + R(u, At) 



(2.5) 
(2.6) 



\R(u,At)\ < const • (At) 2 
Substituting this into (2.5) and taking the limit as At — > gives 

g(u) djj, (u) = 



(2.7) 

Our next step is to derive (2.6) for small At > and compute g(u) explicitly. Let 
X(t) = X and W(t) = u. Due to our assumption (1.4) the probability that the domain 

Di = {(x,v): < (X - x) ■ sgnv < VAt, v min < \v\ < w max } 

contains more than one atom is 0((At) 2 ), hence we can ignore it and assume that D\ 
contains at most one atom. Note that Di is much larger than D defined by (2.2). Now 
it is easy to see that the piston experiences at most one collision during the interval 
(t, t + At). Precisely, a collision occurs if and only if the domain D = D + U D~ contains 
an atom. When it does, and the atom has velocity v at time t, then the piston's velocity 
at time t + At is computed by (1.1): 



W(t + At) 



2m 



u 



M + m 

Therefore, the conditional expectation E(u, At) is 

2m 



(v — u) 



E(u, At) = u + 



M + m 



' vpdxdv — u pdxdv 
d Jd 



0{{Atf) 



(2.8) 



(2.9) 



(the last term comes from our assumption that D contains at most one atom). While 
the second integral in (2.9) is already computed in (2.4), the calculation of the first one 
is similar. For every k > we have 



D 



v h p(x : v)dx = At 



Thus, 



E(u,At) =u + 



At 


'P a V 

--' "mm 


At 




At [Q k+1 - Q 


u + 


2m 


M + m 



(y — u) p-(v) dv + 



v k {u — v) p+(v) dv 



(2.10) 



Q 2 - 2Q lU + Q u 2 ] At + 0{(At) 2 ) 
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This gives us 

9( u ) = tt-, — \Qz - 2 Q± U + Qou 2 

M + m 1 

We now apply the key identity (2.7) and arrive at 

Q2 - 2Q 1 (u) + Q (u 2 ) =0 

where (•) means averaging with respect to the stationary measure /j,q (we note that Qi 
do not depend on u). Thus, the average velocity of the piston is given by 

m _9-m±Q> ( ,„) 

The mechanical equilibrium is defined by the equality of pressures on both sides of 
the piston, i.e. by 

Q 2 = F 2 - F+ = 

and we will restrict ourselves to this case in the rest of the paper. Then (2.11) reduces 
to 

(ho = 9*m M 

which is still different from zero, albeit very small. 

It is well known that at mechanical equilibrium the average value of (W 2 ) is 0(m/M), 
see, e.g. [H, DGL]. But we will estimate it more accurately below. 

Similarly to E(u, At), we define 



E 2 (u, At) = J w 2 V(u,dw;At) 



the conditional expectation of the square velocity W 2 (t + At) given that W(t) = u. The 
stationarity of the distribution /x implies 

J E 2 (u, At) dfj, (u) = J u 2 dfx (u) 

Assume that 

E 2 (u, At) = u 2 + g 2 {u) At + R 2 (u, At) 
with \R 2 (u,At)\ < const • (At) 2 . Then, as in (2.7) 

J g 2 {u) dfi {u) =0 (2.13) 

Now, squaring the equation (2.8) we obtain 

W 2 (t + At) =u 2 + J m (v - u)u + f J m (v - u) 2 
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Then, as in (2.9), 

4m 



E 2 (u,At)=u 2 + — — — J (v — u)up dx dv 

^/ D (,-^«, + ((Af) 



+ 



Computing these integrals as before gives 

4/77, 

g 2 (u) = — (Q 2 u - 2Q lU 2 + Q u 3 ) 

M + m 

+71^ ^ (^3 - 3Q 2 « + 3Q lM 2 - Q u 3 ) 

(M + m) 2 

We included the terms with Q 2 for the sake of clarity, even though we had assumed that 
Q 2 = 0. In the subsequent expressions we remove all Q 2 s. The equation (2.13) gives 

Til 

2Qi(W 2 ) - Q (W 3 ) = —— (Q s + SQ.iW 2 ) - Q (W 3 )) 

M + m 

For brevity, we denote e = m/M. By making one step further and analyzing E^(u, At) = 
J w 3 V(u,dw; At) in the same manner one can see that (W 3 ) = 0(e 2 ), we omit details. 
Hence we arrive at 

r2\ , srs( 2\ 



{w) = Wi £ + oi ] (2 ' 14) 



Substituting this into (2.12) yields 



(W) = ^e + 0(e 2 ) (2.15) 



3 Thermo dynamical analysis 

In this section we express our rigorous results in statistical mechanical terms, such as 
temperature, density, pressure, and heat transfer. 

First, we consider the heat flow across the piston, from right to left. Even though 
both gases are infinite, the heat transfer from one to the other is well defined and its rate 
can be computed. The thermodynamical definition of the heat transfer [GF] is 

R + ^_ = (dS-/dt) - (dM-/dt)(W) (3.1) 

where dS-/dt is the rate of change of kinetic energy and dM.-jdt the rate of change of 
momentum on the left hand side. By the conservation of energy and momentum during 
collisions, it is enough to compute the average change of energy and momentum of the 
piston when it collides with an atom on the left hand side. 
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According to (1.1), the momentum of the piston at one collision changes by 

M + m 

The conditional average change of momentum during the interval (t,t + At), given that 
W(t) = u, is 



1m M r 

E U (SM M ) = 77 / (v-u)dxdv + 0((At) 2 ) 

M + m Jd- 

2 "' } ' (F 2 - - 2F±u + FqU 2 ) At + 0((At) 2 ) 



M + m 

where only collisions with atoms on the left are taken into account. Negating and aver- 
aging with respect to the stationary measure fiQ gives the average momentum transfer to 
the left gas: 

(dM^/dt) = -2m[F~ - 2F~(W) + F (W 2 )} + 0(m 2 /M) (3.2) 

By squaring (1.1) we find the change of the kinetic energy of the piston at one collision: 

5S M = 2mev 2 + 2mVv - 2mV 2 + ■■■ 

where only essential terms are shown (the rest will eventually end up in 0(m 3 /M 2 ) in 
the expression (3.3) below, so we omit them). The conditional average change of kinetic 
energy during the interval (t, t + At), given that W(t) = u, is (again only collisions with 
atoms on the left are taken into account) 

E U {5S M ) = 2m(eF^ - eF^u + F^u - 2F{u 2 + F^u 3 ) At + ■ ■ ■ + 0{{At) 2 ) 

Negating and averaging with respect to the stationary measure /i gives the average 
kinetic energy transfer to the left gas: 

(dS^/dt) = -2m[eFz + F^(W) - 2Ff {W 2 }} + 0(m 3 /M 2 ) (3.3) 

Combining (3.2) and (3.3) we obtain the expression for the heat transfer 

= -2m[eF 3 - - 2F{ (W 2 )} + 0(e 2 ) (3.4) 

Substituting (2.14) gives 

R^- = ^(KF 3 - - FfF 3 + ) + 0(e 2 ) (3.5) 

Now, even though the thermodynamical temperatures of the gases are not defined, 
unless their velocity distributions are Maxwellian, we can define the "effective tempera- 
tures" via the second moment of the velocity distributions: 

rp ,-i fv 2 p±dv x F 2 ± 

T ± = k B m — — = k B m-=± 

Jp±dv F 
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where k B is Boltzmann's constant. We recall that the pressures of the gases are given 
by P± = 2mF 2 ± and their spatial densities by n± = 2F ± , hence the classical law P± = 
n±k B T± holds. 

The piston's velocity distribution /i is, generally, not Maxwellian, but its effective 
temperature can be computed similarly: 

T M = k~ B l M{{W 2 ) - (WO 2 ) ~ k^m- Q> 



2Qi 

where we used (2.14) and (2.15). Hence T M is of the same order of magnitude as T ± , but 
for generic distributions p± there is no simple relation between the exact values of these 
temperatures. 

However, we can derive some interesting relations assuming that the distributions p± 
are Maxwellian. Since this would violate our assumptions (1.4), our further analysis is 
heuristic. Assume that the densities p± satisfy 

p±(x,v) = n±f±(v) 

where n± are (constant) spatial densities and 

1 ( v 2 \ 

f±{v) = i = exp 



2nai V 2 4 

the Maxwellian distributions with temperatures T± = k^ma 2 .. As before, we assume 
that the pressures are equal, i.e. n_T_ = n + T + . It is then easy to compute 

Qo = t, ( n - ~ n +) 

Qi = -}= (n-*- + n + a + ) = ^S= (n^J 2 + n+T] 12 ) 
y2n v27rm v 7 



Qz = \ ( n - (7 - ~ n + a l) = 7^ (™- T - - n + T + ) = 



2k 



Qs = \ - (n_oi + n + a%) = ^ (n^J 2 + n^ 2 ) 
V 7T Virm s v 7 

Now the effective temperature of the piston is 



T - mQ3 - n - T ~ 2 + n+T + /2 - FTF (16) 



2k B Q\ n^T'J 2 + n + Tl /2 
where we used the equality of the pressures. The average velocity of the piston is 

^/27ck B m (n_ - n + )(T_T + ) 1 / 2 



\j2ixk B m 
AM 



T+ - ^/T_ ) + 0{e 2 ) (3.7) 
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Lastly, the heat transfer across the piston (from right to left) given by (3.5) is now 




(3.8) 



Hence, the heat flow is proportional to the temperature gradient, T + — T_, and the heat 
conductivity is thus 



The formulas (3.6)-(3.9) were obtained earlier in [GP, GF] by means of kinetic theory. 

Let us examine closely what happens when the pressures are equal, P — P + , but 
the temperatures are different. Without loss of generality, assume that T_ < T + , i.e. 
the gas on the left is cooler but denser and the gas on the right is hotter but sparser. 
Then (W) > 0, i.e. the piston slowly drifts in the direction of the hotter gas. At the 
same time > 0, i.e. the heat slowly flows from the hotter gas to the cooler gas. 

This demonstrates that the mechanical equilibrium is not a stable state. If our gases 
were finite, there would be a continuing evolution toward thermal equilibrium, in which 
the temperatures become equal, too. Of course, this conclusion makes no sense in the 
idealized model of infinite gases. 

We note that the evolution of the system is somewhat counterintuitive. When the 
piston moves to the right, i.e. V(t) > 0, then by (1.2) the atoms on the right bounce off 
it with a higher speed and so gain energy, while the atoms on the left collide with the 
piston and slow down, hence lose some energy. When the piston moves to the left, i.e. 
V(t) < 0, it is vice versa. Since, on the average, the heat flows from right to left, one 
may conclude that the piston's movements to the left dominate. On the other hand, the 
piston slowly drifts to the right, so its displacements in that direction actually dominate. 

To explain this "paradox", we first recall that the average speed of the piston is of 
order (W 2 ) 1 ^ 2 ~ e 1//2 , which is much larger than the speed of the drift (W) ~ e. Hence 
the piston jiggles back and forth much faster than it drifts in one direction. The heat 
flow is due to "jiggling" rather than "drifting". Indeed, when the piston jiggles to the 
right, the slow atoms on the left have less chance to collide with the piston, since the 
relative velocity is small. Most of the collisions between the slow atoms on the left and 
the piston occur when the latter jiggles to the left. This is why the cooler atoms speed 
up and gain energy, on the average. On the other hand, the hotter atoms on the right 
are less sensitive to the variations in the piston's velocity, and the reason why they cool 
down is different, as we explain next. 

The piston's vibrations are not spatially symmetric. During excursions to the right, 
relatively few collisions with atoms on the right occur, since they are more energetic and 
able to quickly reverse the velocity of the piston. On the contrary, when the piston drifts 




(3.9) 
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to the left, the atoms on the left are weak and it takes them longer to turn the piston 
back to the right. During these intervals, the fast atoms on the right continue hitting 
the piston and lose energy. This explains why the hotter atoms mostly slow down. We 
must admit though that the precise mechanism of the heat flow across the piston in our 
model remains unclear. Some physicists call it a "conspiracy" between the microscopic 
vibrations of the piston and the incoming atoms of the gases [GF, GP]. 



4 Justification of Markov approximation 

Here we estimate the difference between the true velocity of the piston V(t) and our 
Markov approximation W(t). 

Recall that (W 2 ) ~ e = m/M, i.e. the typical speed of the piston (in the Markov 
approximation) is of order y/e. The following theorem, which is essentially proved in 
[CLS] (see remarks below), estimates large deviations of the true velocity V(t): 

Theorem 4.1 ([CLS]) Let \V(0)\ < e 1/2 . There are constants C, d > such that for 
every T > we have 

P { sup \V{t)\ > Ce 1/2 \iie- 1 ) < T 2 e~ dlnlnE 

\0<t<T J 



This gives a good bound on large deviations during time intervals of length T ~ e~ A 
for any large but fixed constant A > 0: it shows that the true velocity V(t) of the piston 
remains 0{e x l 2 lne -1 ) with "overwhelming" probability. 

In our construction of the Markov approximation W(t) we can choose any small but 
positive cutoff value V > 0, and then assume that e is so small that C e 1 ^ 2 hie -1 < V. 
Then we arrive at 

Corollary 4.2 For all sufficiently small e and allT > we have 
P(v(t) = W(t) Vt G (0,T)) > 1 - T 2 e- dhlln£ 

Therefore, during time intervals of length T ~ e~ A for any fixed constant A > 0, 
the true velocity V(t) of the piston coincides with its Markov approximation W(t) with 
"overwhelming" probability. If V(t) and W(t) differ, though, then we have an obvious 
bound: \V(t)\ < tWx- So we obtain analogues of our main formulas (2.15) and (2.14): 

Corollary 4.3 Lets > be small and \ V(0)\ < s 1 ' 2 . Then for any fixed A > and all 
< t < T = s~ A we have 

E(V(t)) = ^e + 0(e 2 ) 

and 



E(V 2 (t)) = ^-e + 0(e 2 ) 



where E(-) is the mean value. 
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It follows from the results of [CLS] that whenever the true velocity V(t) exceeds 
C e 1 ! 2 lne -1 , then with "overwhelming" probability it will be driven back toward zero by 
further collisions with the gas atoms. Hence, the estimates of the last corollary should 
remain true for all t > 0, but we do not pursue such a goal, because the drift of the piston 
becomes evident on the time scale of order e~ A already for A > 2, as we show next, so 
the estimates in Corollary 4.3 are sufficient for practical purposes. 

Indeed, let X(0) = and consider the random position of the piston X(t) at time t. 
Also, let Y(t) = Jq W(s) ds be the corresponding function in the Markov approximation 
model. We obviously have 



(Y(t))=t(W) = ^et + 0(e 2 t) (4.1) 

We also need to estimate the variance of Y(t), which is a little harder. We claim that 

Var[F(t)] < Dt (4.2) 

with some constant D > independent of e. Below we outline the argument, suppressing 
some technical details that are based on the estimates obtained in [CLS]. 
It is standard that 

Var[y(t)] = /* f Cov(W(s),W(u))dsdu 
Jo Jo 

By the stationarity of the process W(t), the covariance here equals 

Cov(W(s),W(u)) = p s „ M Var(W/) 

where p s _ u is the correlation between W(s) and W(u), which only depends on \s — u\. 
Also recall that 

Var(W0 = (W 2 ) - (W) 2 = e + 0(e 2 ) 

2Qi 

Therefore, 

Var[y(t)l < const -te [ p s ds (4.3) 

Jo 

It remains to estimate the decay of the velocity autocorrelation function p s . 

Suppose during the interval (u, u + s) the piston experiences k collisions with gas 
atoms with velocities vi,...,Vk numbered in the order in which the collisions occur. 
Then (1.1) implies 

k 

W(u + s) = (1 - a) k W(u) + a ^(1 - a) k ' j ■ Vj (4.4) 

j=i 

where a = 2m/ (M + m) ~ 2e. This equation can be easily verified by induction on k. 
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Now, the number of collisions k grows as const-s, on the average, in fact 

P(\k- Cl s\ > c 2V ^lns) < s 2 £-d^ (4.5) 

with some positive constants Ci,C2, see [CLS]. For s < e~ A this gives a good bound on 
large deviations, so we can safely assume that k > C\sj2. Next, the velocities i>j are 
almost independent of W{u). In fact, there is a dependence, since the time of collision 
between a particular atom and the piston depends on the velocity of the latter, but the 
correlation between W(u) and velocities Vj for every fixed j is negative. This follows 
from a simple observation: if W(u) > 0, the piston is more likely to collide with atoms 
on the right, for which v j < 0, and vice versa. Therefore, the correlation p s is determined 
by the first term in (4.4), hence 

\p a \ < (1 - «) C1S/2 

This implies 

rt roo 9 I 

/ p s ds< / (1 -a) cis/2 ds = — r~ — 

Jo Jo Ci m(l — a) C\E 

and combining with (4.3) proves (4.2). 

The estimates (4.1) and (4.2) show that the piston's displacement due to the drift 
exceeds its typical random fluctuations at times t = e A for all A > 2. On such time scales 
the drift is "physically evident" and it is possible to verify it experimentally. Numerical 
experiments of this sort were done by [GF]. 

We now comment on the proofs of Theorem 4.1 and the estimate (4.5). The corre- 
sponding results were proven in [CLS] in a slightly different context. That paper dealt 
with a heavy piston in a cubical container of side L filled by an ideal gas. That model, 
too, reduced by a trivial projection onto the axis perpendicular to the piston surface to 
a one-dimensional gas on an interval [-L/2, L/2]. The velocity cutoff assumptions in 
[CLS] were identical to our (1.4). The piston had mass M = aL 2 (a > was a constant) 
and the spatial density of the one-dimensional gases was proportional to M — aL 2 , but a 
simple rescaling of time and space by M reduces the spatial density to values 0{1) and 
makes it independent of M. Hence, we arrive exactly at the same initial conditions as 
in the present paper, except the gases in [CLS] were finite - they evolved on an interval 
[— aL 3 /2, aL 3 /2] with the piston initially placed at X(0) = 0. 

Due to the velocity cutoff (1.4) the gas atoms in [CLS] cannot interact with the piston 
more than once unless they travel across the half interval [— aL 3 /2, aL 3 /2], reflect at an 
endpoint x = ±aL 3 /2 and then travel back to the piston in the middle. This takes 
time T = 0(L 3 ), as it was proven in [CLS]. Thus, during the initial interval (0,T) with 
T = 0(M 3 / 2 ) the evolution of the piston studied in [CLS] is identical to ours. Hence, 
all the results and estimates derived in Section 3 of [CLS] directly apply to our present 
context. Furthermore, with some minor obvious changes those results can be extended 
to arbitrary time intervals as we did in Theorem 4.1 and (4.5). 

Lastly, we have assumed the velocity cutoff (1.4) here and in the paper [CLS] mostly 
for convenience. It might be true that the upper bound on velocities is quite essential 
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in [CLS], since arbitrarily fast atoms could bounce back and forth between the piston 
and the wall many times during short time intervals thus making significant impact on 
the piston, as it was remarked in Section 5 of [CLS]. But in the present work, if we 
had arbitrarily fast particles, they would hit the piston once and get away never to come 
back, so their impact would be quite limited. 

As for the lower bound on velocities, it seems to be less essential in both models. It 
was conjectured in Section 5 of [CLS] that the results of that paper could be extended 
to velocity distributions without lower bounds (i.e. with t> min = 0). If so, then the 
estimate in Theorem 4.1 implies that the typical velocity of the piston remains of order 
s 1 / 2 \ris~ 1 . As long as this is true, only atoms with smaller velocities can experience two 
or more collisions with the piston, thus violating the Markovian character of the piston's 
velocity process V(t). It is customary to redefine the dynamics so that the collisions 
with slow atoms are ignored altogether, and then the resulting process W(t) will be 
Markovian. The difference \V(t) — W(t) \ should be estimated separately and presumably 
is negligible. Such a strategy was successfully implemented in earlier works [H, DGL] in 
somewhat different context, and it is likely to go through in our case as well. This is the 
subject of our future work. 
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